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We analyze the atmospheric neutrino experiments of SuperKamiokande (535 days) using 
CTs ■ the three-flavor neutrino framework with the mass hierarchy rrt\ ~ m% <C 7713. We study 

the event ratios of sub-GeV, multi-GeV and upward through-going muons zenith angle 
distributions. Taking into account the atmospheric and terrestrial experiments, we obtain 
the allowed regions of mass and mixing parameters Am\ 2 , sin 2 2#i2, Z\m 2 3 , #13 and $23- 
The mass parameter AmW is restricted to 0.01 eV 2 -0. 0002 eV 2 . For mixing parameters, 
there is no difference between the large Q\i angle solution and the small one, and #13 < 
13°, 28° < 6*23 < 37° and 53° < <9 23 < 62° for Am\ z = 0.01 eV 2 , 9 13 < 23° and 29° < 6» 23 < 
61° for Am\ z = 0.002 eV 2 , 4° < 13 < 22° and 38° < 6> 23 < 54° for Am\ z = 0.0002 eV 2 . 
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O . §1. Introduction 

G\ 

SuperKamiokande experiment HP has observed a definite atmospheric neutrino anomaly 
and confirmed a hypothesis of neutrino oscillation. In their two-flavor mixing analyses of the 
sub-GeV and multi-GeV zenith angle distribution, it has been obtained that the <-» v T 
oscillation is preferred over the «-> v e oscillation and the range of mass parameter Am 2 is 
from 3 x lCr 4 eV 2 to 0.8 x lO^eV 2 .^ 

Observing the neutrino oscillations is one of the most important experiments suggesting a 
new physics beyond the standard model. Constructing a model beyond the standard model, 
it is necessary to confirm the mixing and mass parameters of three neutrinos in the low 
energy region although there is no restriction on number of neutrinos in very high energy 
region, i- 1 The three- flavor neutrinos scenario!!* is found to be consistent with above the large 
Uu <-» u T oscillation and the solar neutrino anomaly^ and terrestrial neutrino experiments. 

K2K, MINOS and CERN-Gran Sasso long-baseline experiments!^ are in preparation 
to confirm the <->• v e and <-> v T oscillation precisely. Then it becomes important 
to examine the mass and mixing parameters in the three-flavor neutrinos framework&'& 
precisely. 

In this paper, we analyze the SuperKamiokande atmospheric neutrino experiment (535 
days) of the sub-GeV, multi-GeV and upward through-going muons zenith angle distribu- 
tions in the three-flavor neutrino framework with a hierarchy of neutrino masses rri\ ~ 
m 2 7713. In a previous work,© we have analyzed the zenith angle distributions in the 
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SuperKamiokande atmospheric neutrino experiment(325.8 days). HIP In that work, we ana- 
lyzed the double ratios of the sub-GeV and multi-GeV zenith angle distributions neglecting 
the matter effects of the Earth and the difference between incident neutrinos and detected 
charged leptons. 

In the three- flavor neutrino framework with the hierarchy mi ~ m 2 <C m 3 , there are 5 
parameters (9 12 , #13, 623, Am\ 2 = m\ — m 2 , Am 23 = m 2 — m\ concerned with the neutrino 
oscillation. For solar neutrino deficit, the MSW solutions' 1 '©'© taking into account the 
matter effects predicts the large angle solution sin 2 29\ 2 = 0.6-0.9, Am\ 2 = 4 x 10~ 6 - 
7 x 10" 5 eV 2 and small angle solution sin 2 29 12 = 0.003-0.01, Am\ 2 = 3 x 10~ 6 -1.2 x 10~ 5 eV 2 
for 9 1 3 = 0°-20°, and these large and small angle solution are merged for # 13 = 25°-50°. 
The vacuum solution for solar neutrinos is obtained as Am\ 2 ~ lO _lo eV 2 .0' ) In terrestrial 
neutrino experiments, Fogli et a/.© have thoroughly analyzed and obtained the allowed 
regions in tan 2 # 13 -tan 2 9 23 plane for various values of the mass parameter Am 23 . We also 
analyzed the terrestrial neutrino experiments, and obtained results similar to those of Fosrli 
et al. The allowed region in tan 2 #13 — tan 2 9 23 plane was more restricted for higher values of 
Am 23 , and the allowed region spread out for lower values of it. In this paper, we combine 
the results obtained from the atmospheric neutrino experiments of SuperKamiokande and 
the terrestrial experiments to obtain the allowed regions for mixing and mass parameters in 
three-flavor neutrinos. 



§2. Neutrino oscillation in three- flavor neutrinos 

The unitary matrix U representing the neutrino mixing is defined as 

3 

v i a = 2 u i*P v P> l <* = e ' ^ r ' ( X ) 

/3=1 

where the states v\ a and up are the flavor and mass eigenstate of neutrinos, respectively. This 
mixing matrix corresponds to the CKM matrix Uq KM for quarks sector. We parameterize 
the unitary matrix neglecting the CP violation phases as 

U = exp (i9 23 \ 7 ) exp (i6> 13 A 5 ) exp (i9 12 \ 2 ) 
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where the \ are the Gell-Mann matrices. 

The probabilities for transitions v\ a — > v\ & are written as 

P(yi a -> %) = I < viM^M > I 2 = +pll-+v lB Si2 
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where is the term representing the neutrino oscillation; 



/\ 177,?. 

= sin 2 1.27— ^L. (4) 
E 

Here /^m.^- = |m 2 — m 2 |, E and L are measured in units eV 2 , GeV and km, respectively. 

The values of neutrino masses are not known precisely, but two mass parameters are 
necessary to account for the solar neutrino experiments and the atmospheric neutrino ex- 
periments. In the former, mass parameter Am 2 is obtained as 10 _4 -10 -5 eV or 
~ l()- lo eV 2 ,0) and in the latter, Am 2 is obtained as 10 -1 - l(T 3 eV 2 .0)'@) Then it seems 
that the most reasonable mass hierarchy is that in which the lower two neutrino masses of 
three neutrinos are very close and the third is rather far from them. Thus we assume that 
three neutrino masses have a mass hierarchy obeying 

mi Pd m 2 <C 7713. (5) 

With the hierarchy Eq. (||), Am 2 2 <C Am^ ~ Am 2 3 , Eq. (3) for the transition probabilities 
P{ u i a ~~ * v i&) can be rewritten as 

Piyu - n.) = 1 - 2(1 - 2Ul 3 - uf al - uf a2 + UfJS 12 - 4f/ 2 3 (1 - uf a 3 )S 23 , 

+AUl 3 U? g3 S 23 . (6) 

In the atmospheric neutrino experiments, especially in the sub-GeV ones, the matters 
of the Earth have the important effects. Matter effects are induced by the quantity A = 
2\f2EGpN e . In the Earth, N e ~ 3mol/cm 3 , thus ylsin2# 13 << 1 and we can approximate 
the mixing matrix U as© 

U M = exp(iA 7 23 ) exp(a 5 13 ) exp(zA 2 (7) 

where 

sin2^| = ^sin2^ 12 , 

Am™ 2 = mf 2 - mf 2 , E = m\ + m 2 , 
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\{{E + Acos 2 6 13 ) =f sJ(A cos 2 9 13 - Am\ 2 cos 2fl 12 ) 2 + (Am 2 12 sin 2# 12 ) 2 } 



tti^ 2 = m| + A sin 2 6»i 3 . 
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This expression corresponds to the neutrino case, and A enters in this expression with 
opposite sign for the anti-neutrino case. If we approximate the density of the Earth to be 
uniform, we can take the expressions for the transition probability with matter effects as 




<M 

23 J 



(8) 



where 



Sff = sin 2 1.27 




(9) 



3. Numerical analyses of neutrino oscillations 



3.1. Atmospheric neutrinos 

Evidence for an anomaly in atmospheric neutrino experiments was pointed out in reports 
of Kamiokande collaboration and 1MB collaboration usingwater-Cherencov exper- 
iments. More recently, reports of SuperKamiokande collaboration E-^'ElI' have provided more 
precise results on the anomaly in atmospheric neutrinos. In our previous papers, 0) we have 
analyzed the double ratios i?(/i/e) = -R expt (/i/e)/i?Mc(/ i / e ) of atmospheric neutrinos. In this 
work, we analyze the event ratio N Exp (l a ) / N M c(l a ) , where l a represents /i and e. 

The zenith angle distributions of sub-GeV, multi-GeV events and upward through-going 
fx fluxes are taken from the SuperKamiokande 535 days experiments. The data are tabu- 
lated in Table I. These values are taken from the experimental event data and Monte-Carlo 
simulations which are given graphically in Ref. 1). /x-like events include the fully contained 
and partially contained events. Errors represent statistical ones only. In the above data, 
sub-GeV experiments detect the visible-energy less than 1.33GeV. 

The zenith angle 9 dependent events dN Exp (l a ) / d cos 9 and dN MC (l a ) / d cos 9 are defined 

as 



where the summation Y^a are taken over and v e . In these expressions, processes of 
andz/ e are contained. e a (E a ) is the detection efficiency of the detector for a-type charged 
lepton with energy E a . a a (E u , E a ,ip) is the differential cross section of scattering l a with 
energy E u by incident v a with energy E a , where angle ip is the scattering angle relative to 
the direction of incident v a . Fp(E v ,9) is the incident vp flux with energy E v at the zenith 
angle 9. P M {y a — > z/3) is the transition probability with the matter effects expressed in 



dN Exp (l a ) 



= J2 J e a {E a )a a {E v ,E a ^)F p {E v1 9 - ^)P M {u B - v a )dE a dE v d cos ijj, (10a) 



dcos9 




= j e a (E a )a a {E v ,E a ^)F a (E v1 9 -^)dE a dE v dcos^i (10b) 
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Table I. e-like, /x-like atmospheric neutrino data and upward through-going muons fluxes of Su- 
perKamiokande experiments (535days data). These values are the ratios of experimental data 
and Monte-Carlo simulations which are obtained from the graphs in Ref. 1). ii-like events 
include fully contained and partially contained events, and errors represent statistical ones 
only. 



Sub-GeV data 


cos 9 


range 


e-like event ratio N Exp /N MC fJr 


-like event ratio N Exp /N MC 


-1.0 • 


- -0.6 


1.37 ±0.09 


0.56 ±0.04 


-0.6 


- -0.2 


1.12 ±0.08 


0.71 ±0.05 


-0.2 


- 0.2 


1.18 ±0.08 


0.74 ±0.05 


0.2 • 


- 0.6 


1.05 ±0.08 


0.86 ±0.06 


0.6 


- 1.0 


1.15 ±0.08 


0.82 ±0.05 


Multi-GeV data 


COS 6 


range 


e-like event ratio N Exp /N MC fJr 


dike event ratio N Exp /N MC 


-1.0 


- -0.6 


1.35 ±0.21 


0.56 ±0.07 


-0.6 


- -0.2 


1.10 ±0.16 


0.57 ±0.07 


-0.2 


- 0.2 


1.13 ±0.15 


0.79 ±0.07 


0.2 


- 0.6 


1.42 ±0.18 


1.02 ±0.09 


0.6 


- 1.0 


1.18 ±0.20 


1.03 ±0.10 



Upward through-going /x data 
cos 6 range \i flux ratio N Exp /N MC 



-1.0 - 


-0.9 


0.77 ±0.14 


-0.9 - 


-0.8 


0.79 ±0.12 


-0.8 - 


-0.7 


0.58 ±0.11 


-0.7 - 


-0.6 


0.96 ±0.13 


-0.6 - 


-0.5 


0.73 ±0.10 


-0.5 - 


-0.4 


0.80 ±0.10 


-0.4 - 


-0.3 


0.74 ±0.10 


-0.3 - 


-0.2 


0.92 ±0.10 


-0.2 - 


-0.1 


1.02 ±0.10 


-0.1 - 


0.0 


1.09 ±0.10 
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Eq. (8) and it depends on the energy E v and the distance L. This distance depends on the 
zenith angle 9 as L = 



(r + h) 2 — r 2 sin 2 9 — r cos 0, where r is the radius of the Earth and 
h is the altitude of the production point of atmospheric neutrinos. 

Informations of Fp(E u , 9) for multi-GeV case are given in Refs. ^),|21|),^3D. Informations 
of Fp(E v ,9) including the geomagnetic effects for sub-GeV case are taken from Ref. |22| ). 
Other informations of e a (E a ) and a a (E u , E a , ip) are given in Ref. [24|). The upward through- 
going muons fluxes are simulated in Ref. |23J). Explicit calculation of Eq. (10) is explained in 
Appendix A. 
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(c) (d) 
Fig. 1. The plots of allowed regions in the tan 2 #i3~tan 2 #23 plane determined by the zenith angle 
distributions of SuperKamiokande 535 days data. These figures correspond to the large #12 
angle solution, Am\ 2 = 3 x 10 _5 eV 2 and sin 2 29 = 0.7. In these figures, broken thick, broken 
thin and dotted curves denote the regions allowed in 99%, 95% and 90% C.L., respectively. 
Figs, (a)-(d) show the plots of sub-GeV experiments. 



Since P M {y a — > up) is the function of Am\ 2l Am 23 , 9u, #13 and #23, the ratio (dNEx P (l a ) 
I 'dcos9) I '(dNyic(la) I 'dcos9) of the zenith angle distributions is a function of Am\ 2) Am^, 612, 
#13, #23 and 9. We analyze the atmospheric neutrino data fixing the parameters Am\ 2 
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Fig. 1. The plots of allowed regions in the tan 2 #i3~tan 2 023 plane. Figs, (e)-(h) the plots of multi- 
GeV experiments. 

and sin 2 2^12, which have been determined from the solar neutrino experiments, i>'lll)'0) 
as Am\ 2 = 3 x lCT 5 eV 2 and sin 2 2#i2 = 0.7, which corresponds to the large angle solution, 
and Am\ 2 = 10 _5 eV 2 and sin 2 2^12 = 0.005, which corresponds to the small angle solution. 
We treat the ratios of the zenith angle distributions of the experimental events and the ones 
of Monte-Cairo simulation, (iV Exp (Z Q )/iV MC (/ Q ,))j, where i represents the region number of 
the bins of zenith angle 9. 
We treat the \ 2 defined as 



X 



E 



{(N Exp (l a )/N MC (l a ))f - (N Exp (Q/N MC (l a ))t^Y 



For the sub-GeV and multi-GeV experiments, the summation on % are over 1 to 5 of zenith 
angle range bins and for the upward through-going fi, 1 to 10. The summation on l a are 
over \i and e for the sub-GeV and multi-GeV experiments. We estimate the \ 2 on the 
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Fig. 1. The plots of allowed regions in the tan 2 #i 3 -tan 2 #23 plane. Figs, (i)-(l) show the plots of 
the combinations of the sub-GeV and multi-GeV experiments. 



various of Z\m| 3 ,0i3 and 6*23. cr st represents the statistical error and a sy systematic one. 
We adopt the values of a sy which are induced from the magnitudes of 10% of the Nuc 
assumed as uncertainty. These values are deduced from the graphs given by the report 
of SuperKamiokande collaboration. 0* In Fig. 1, we showed the contour plots of x 2 i n the 
tan 2 # 13 -tan 2 #23 plane for various Am^ z . We showed the plots of sub-GeV experiment in 
Fig. l(a)-(d) and the plots of multi-GeV one in Fig. l(e)-(h) for the large 8 12 angle solution, 
Am\ 2 = 3 x 10" 5 eV 2 and sin 28 = 0.7. In these figures, broken thick, broken thin and 
dotted curves denote the regions allowed in 99%, 95% and 90% C.L., respectively. The 
combinations of the sub-GeV and multi-GeV experiments are shown in Fig. 1 (i)-(l). Fig. 2 
represents the case of the small #12 angle solution, Am\ 2 = 10 _5 eV 2 and sin 2 28 = 0.005. 

From these plots, we can get the following results for the mixing parameters: 
(1) Am\ 3 is allowed from 10 _2 eV 2 to 2 x 10 _4 eV 2 . The minimum x 2 is obtained as 14.5 
at Aml 3 = 7 x 10 _4 eV 2 , tan 2 13 = 6 x 10~ 2 and tan 2 823 = 2. The minimum x 2 m tne 
restriction 6*13 = 0, which corresponds to the — v T mixing, is 17 at Am^ = 2 x 10 _3 eV 2 . 
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Fig. 2. The plots of allowed regions on tan 2 #i3-tan 2 #23 plane determined by the zenith angle 
distributions of the combination of sub-GeV and multi-GeV fi SuperKamiokande 535 days data. 
These figures correspond to the small #12 angle solution, Am\ 2 = 10 _5 eV 2 and sin 2 29 — 0.005. 
In these figures, broken thick, broken thin and dotted curves denote the regions allowed in 99%, 
95% and 90% C.L., respectively. 



(2) v e — v T mixing is small and — v T mixing is large; tan 2 # 13 < 10 _1 and tan 2 #23 ~ 1. 

(3) There is no significant difference between the large 9yi solution and small one. 

The result that Am^ is allowed from 10 _2 eV 2 to 2 x 10~ 4 eV 2 and the minimum \ 2 m the 
restriction #13 = is obtained at Am\ z = 2 x 10 _3 eV 2 is the same as the result obtained by 
SuperKamiokande collabolation. 0) 

Plots of the upward through-going \i are shown in Fig. 3(a)-(d) and the combinations of 
the sub-GeV, multi-GeV and upward through-going \x flux in Fig. 3(e)-(h) for the large 9 12 
angle solution. The plots for the small 9 12 angle solution is the similar to the Fig. 3. From 
these result, we can say that 

(1) the mixing parameters allowed in the sub-GeV and multi-GeV zenith angle distributions 
can explain the upward through-going fi flux experimental data, 
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(2) in large value of Am\ z ~ 10 -1 -10~ 2 eV 2 , the values near 45° of 6*23 are excluded, because 
the deficit from 1 of the event ratio Ne xp {[i) / Nyic(fi) is not so large (about 0.2) in upward 
through-going fi flux compared with multi-GeV experiments. 
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Fig. 3. The plots of allowed regions in the tan 2 #i3~tan 2 #23 plane determined by the zenith angle 
distributions of the upward through-going fj, SuperKamiokande data for large #12 angle solution. 
There is no difference between the large 612 angle solution and small one. In these figures, the 
broken thick, broken thin and dotted curves denote the regions allowed in 99%, 95% and 90% 
C.L., respectively. Fig. (a)-(d) show the plot of the upward through-going fi experiments. 



3.2. Terrestrial neutrinos 



By many terrestrial experiments, the strong restrictions are imposed to neutrino mixing 
parameters and masses. In terrestrial experiments, there are two types of experiments, short 
baseline and long baseline experiments.^' In this paper, we analyze the short baseline exper- 
iments. Accelerator experiments searching for the appearance of v T in were performed by 
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Fig. 3. The plots of allowed regions in the tan 2 #i3-tan 2 623 plane for the combinations of sub-GeV, 
multi-GeV and upward through-going [i experiments. 



E531, CHORUS and NOMAdH 



P{y» -> v r ) < 2 x \Q- A (90% C.L.) 
L/E ~ 0.02. 



(12) 



Accelerator experiments searching for the —>■ v e and — » v e oscillations were accom- 
plished by E776, KARMEN and LSNdI, 



P(v„ -> Ve) < 1-5 x 10^ (90% C.L.), E776 

L = 1km, E ~ lGeV, 
Pfa -> < 3.1 x 10~ 3 (90% C.L.), KARMEN 

L = 17.5m, E < 40MeV, 



P{Vn 



9 e ) = 3.4+ 2 ;° ± 0.7 x 10~ 3 , LSND 
L = 30m, E ~ 36 - 60MeV. 



(13a) 
(13b) 
(13c) 



11 



Accelerator experiment searching for the disappearance of was carried out by CDHSW 
experiment and the nuclear power reactor experiments searching for the disappearance of v e 
were carried out by the BUGEY experiment and CHOOZ experiment. & 

z/ M — disappearance experiment, CDHSW (14a) 
0.26eV 2 < Am 2 < 90eV 2 are excluded for maximal mixing, 

L = 130, 885m, E ~ 3GeV. 
1 - P{p e -> v e ) < 10~ 2 (90% C.L.), BUGEY (14b) 

L = 15, 40, 95m, E ~ 1 - 6MeV. 
Pipe -> Ve) = 0.98 ±0.04 ±0.09, CHOOZ (14c) 

L = 1km, E ~ 3MeV. 

We show the contour plots of allowed regions in the tan 2 # 13 -tan 2 823 plane determined 
by the probability P expressed in Eq. (6) satisfying above experimental data Eqs. (12), ( |T3|) 
(except LSND data (13c)) and ( |i~4|) in Fig. 4. Allowed regions are left- and right-hand sides 
surrounded by curves. The curves represent the boundary of 90 % C.L. of P. We fixed the 
parameters Am\ 2 and Qyi as Am\ 2 = 10 _5 eV 2 and sin 2 2di2 = 0.8, and the various values 
of parameter Am 23 between O.leV 2 and 0.0001 eV 2 . Although we fix the parameters Am 2 2 
and #12 as Am\ 2 = 10~ 5 eV 2 and sin 2 29i 2 = 0.01, the results are not changed, because the 
probability P of a terrestrial neutrino is insensitive to the v e — mixing parameters Am\ 2 
and sin 2 26 , i 2 . The dotted lines show the allowed regions restricted by LSND data Eq. (13c). 
The allowed regions determined by LSND data are very restricted, then hereafter we do not 
consider it. 

3.3. Allowed regions of mixing parameters 

Here, We discuss the allowed regions for mixing parameters Am\ 2 , sin 2 20 12 , Am 23 , # 13 
and 823 satisfying atmospheric neutrino and terrestrial neutrino experiments. Am\ 2 and 
sin 2 2# 12 have been determined in analyses considering the MSW effects of the solar neutrino 
experiments: B*' 0)> 0) 

(Am 2 12 , sin 2 26>i 2 ) 

(4 x 10~ 6 - 7 x 10~ 5 eV 2 , 0.6 - 0.9), large angle solution . 
(3 x 10~ 6 - 1.2 x 10~ 5 eV 2 , 0.003 - 0.01). small angle solution 
for 13 = 0° - 20° 

Observing the allowed regions obtained in terrestrial neutrino experiments (Fig. 4) and 
the allowed regions obtained in atmospheric neutrino experiments including sub-GeV and 
multi-GeV data (Fig. 1 and 2), we obtain the allowed regions consistent with these experi- 
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Fig. 4. The plots of allowed regions in the tan 2 ^13-tan 2 #23 plane determined by P of terrestrial 
appearance experiments of — ► u T , — ► v e , — ► u e , and the disappearance experiments 



of v.. 



and v e 



The allowed regions are left- and right-hand sides surrounded by 



curves. The curves represent the boundary of 90 % C.L. of P. Am\ 2 and sin 2 26\2 are fixed as 
l(T 5 eV 2 and 0.8, respectively. Am\ 3 is fixed to 0.1eV 2 ((a)), 0.01eV 2 ((b)), 0.001eV 2 ((c)) and 
0.0001eV 2 ((d)). Dotted lines show the allowed regions determined by LSND data. 



ments. The allowed regions of 13 and #23 for the large 9 12 angle solution are as follows: 

no allowed region, 
9 13 < 13°, 28° < 23 < 37°, 53° < 23 < 62°), 
(013 < 9°, 28° < 023 < 62°), 
(0 13 < 23°, 29° < 023 < 61°), 
(0 13 < 20°, 29° < 23 < 61°), 
(013 < 24°, 35° < 023 < 55°), 
(4° < 13 < 22°, 38° < 23 < 54°), 
no allowed region. 



for Am\ 3 


= O.leV 2 


for Am 2 23 


= O.OleV 2 


for Am 2 23 


= 0.005eV 2 


for Am\ 3 


= 0.002eV 2 


for Am 23 


= O.OOleV 2 


for Am\z 


= 0.0005eV 2 


for Am\ 3 


= 0.0002eV 2 


for Am\ 3 


= O.OOOleV 2 



(16) 
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For the small 6u angle solution, similar solution is obtained. We cannot see the significant 
difference between the large 6u angle solution and small one from these results. 

If we adopt the 90% C.L. of x 2 for the zenith angle dependence of atmospheric neutrinos, 
the range of mass parameter Am\ 3 is restricted to as follows; 

Am\ 3 = O.OleV 2 ~ 0.0002eV 2 . (17a) 

The value of Am\ 3 at the minimum of x 2 under the constraint of #13 = is obtained near 
tan 2 623 = as 

Z\m 2 3 = 2 x l(T 3 eV 2 with X 2 = 17. (17b) 

Under the no constraint of #13, the value of A\\ at the minimum of x 2 is obtained at 
tan 2 613 = 6 x 1CT 2 and tan 2 623 = 2 as 

Am 2 23 = 7 x l(T 4 eV 2 with X 2 = 14.5. (17c) 

These results (17a) and (17b) are the same as those of SuperKamiokande (Ref. |I|) and similar 
to those of a three-flavor neutrinos analysis (the last paper in Ref. |9|)). 

The zenith angle distribution of upward through-going fi flux can be explained by the 
mixing parameter ranges (16) restricted by the zenith angle distributions of sub-GeV and 
multi-GeV experiments and terrestrial experiments except for the values near 45° of 823 in 
large value of Am\ 3 (~ 10 _1 — 10~ 2 eV 2 ). This exception is because the deficit from 1 of the 
ratio -/V Ex p(/^)/A^ M c(/i) is not so large (about 0.2) in upward through-going /z flux compared 
with multi-GeV events. 

We comment on the present results and ones obtained in our previous analysis. 0) In 
previous study, we analyzed the double ratios of the zenith angle distributions of atmospheric 
neutrinos neglecting the matter effects of the Earth and the smearing effects, and obtained 
the results that the mass of Am\ 3 is rather large and the large 612 angle solution is favored 
than the small one. Analyzing the event ratio in present study, the upper limit of the mass 
parameter Am\ 3 is limited to low value. The difference between large #12 angle solution and 
small one is disappeared by containing the matter effects of the Earth. 



§4. Conclusion 

We analyzed the atmospheric neutrino experimental data of SuperKamiokande % under 
the three-flavor neutrino framework with the mass hierarchy mi « m 2 <C m 3 and obtained 
the allowed regions of the parameters Am\ 2) sin 2 2#i 2 , Am\ 3 , 6*13 and 623. We studied 
the event ratios of the sub-GeV, multi-GeV and upward through-going muons zenith angle 
distributions. From these atmospheric experiments and terrestrial ones, the allowed regions 
of mass parameters Am\ 3 are restricted as 0.01eV 2 -0.0002eV 2 , and the value of Am\ 3 at the 
minimum x 2 under the restriction #13 = is 2 x 10~ 3 eV 2 . For mixing parameters, there is 
no difference between the large 612 angle solution and the small 6*12 one. Obtained results 
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are 9 13 < 13°, 28° < 9 23 < 37° and 53° < 9 23 < 62° for Am\ 3 = O.OleV 2 , 9 13 < 23° and 
29° < 9 23 < 61° for Am 2 3 = 0.002eV 2 , 9 13 < 24° and 35° < 6 23 < 55° for Am 2 3 = 0.0005eV 2 , 
4° < 9 13 < 22° and 38° < 9 23 < 54° for Am\ 3 = 0.0002eV 2 . In our three neutrino analysis 
for the SuperKamiokande atmospheric neutrino experiments, the large — v T mixing is 
favored and the small mixing v e — mixing is allowed at the same time. This result will be 
significant in the long baseline experiments. 



Appendix A 

Explicit calculation of Eqs.(lOa) and (10b) 

In multi-GeV case, the information of the Super-Kamiokande detection efficiency e a (E a ) 
is not known for us, but the detected charged lepton event number f a (E a ) is given in the 
paper of SuperKamiokande collaboration.© The event number f a (E a ) is expressed as 

f a (E a ) = J e a {E a )a a {E v ,E a ,il))F a (E v ,9 - if))dE„d cos ipd cos 9. (A-l) 

We define the effective zenith angle dependent neutrino flux f a {E a , 9) as 

fa(E a ,9) = f a {E a ) / a( f;f F flV (A-2) 
J_ x a cos 9F a (E a , 9) 

where we use the zenith angle distribution given in Ref. [H]) for F a (E a ,9). Furthermore 
we approximate the if) dependence of differential cross section <j a by a smearing Gaussian 
function n(if>) as 

M) = 1 f tan^/2) ) 

vvv 0rtan(^ o /2) (1 + cosif))smip P \ tan 2 (^ /2) J ' v ; 



where the resolution if> is taken as 17°. Then, we estimate dN^ ;xp /dcos9 by the following 
integral, 

[n(if))fp(E ,9-if))P M (up^v a )dE p dcosif). (A-4) 
a cos o J 

In sub-GeV (E u ) = j a a (E u , E a , ip)d cos ipdE a is calculated in Ref. ^4|). e a (E a ) 

is nearly 1 in the sub-GeV region, thus we estimate the detected lepton event number with 
zenith angle dependence f a (E u , 9), approximating 9 = if) in F a {E v , 8 — if)), 

f a (E u ,9) = e a (E u ) J a a (E u ,E a ,if))F a (E u ,9)dE (X dcosif) = e a (E u )a a (E u )F a (E„,9), (A-5) 

where information of F a (E„,9) is taken from Ref. ^2]). Then, introducing the smearing 
Gaussian function n(ip) as the multi-GeV case, we express dN-s xp /d cos 9 as 



dcos9 



n(if))ffs(E u ,9- if))P M {v p -> u a )dE u dcosif), (A-6) 
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where the resolution tpo is taken as 60° for sub-GeV case. 

For upward through-going muons fluxes, event number f a (E a ) of through-going muons 
are estimated in Ref. 24). Thus, using this information and the same treatment as multi-GeV 
case, we calculate the zenith angle 6 dependent events of upward through-going muons. 
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